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Calculations are presented for incoherent J/ ip electroproduction from the deuteron 
at JLab energies, including the effects of J/'0-nucleon rescattering in the final state, 
in order to determine the feasibility of measuring the J/V'-nucleon scattering length, 
or the J/V'-nucleon scattering amplitude at lower relative energies than in previ- 
ous measurements. It is shown that for a scattering length of the size predicted by 
existing theoretical calculations, it would not be possible to determine the scatter- 
ing length. However, it may be possible to determine the scattering amplitude at 
significantly lower relative energies than the only previous measurements. 
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I. INTRODUCTION 

With the upcoming 12 GeV upgrade at JLab, electroproduction of the J/ip at JLab on 
a proton or deuteron will be possible. With the mass of the J /if) being 3.097 GeV, the 
threshold virtual photon energy for electroproduction (at small Q"^) on a single nucleon is 
Vthresh — 8.2 GeV in the LAB frame, and is thus accessible with a 12 GeV electron beam. 
Most of the existing data on J/ ijj photo- and electroproduction is at much higher energy. The 
12 GeV upgrade provides the opportunity to measure J/^ production near threshold [2]. In 
addition, measuring electroproduction on the deuteron provides the opportunity to measure 
the J/'j/'-nucleon elastic scattering amplitude at lower energies than previous measurements, 
if the rescattering of the produced J/ip on the spectator nucleon in the deuteron is non- 
negligible. This paper presents calculations of the incoherent J/ip production amplitude from 
the deuteron near threshold (in order to determine if the J/'?/'-nucleon scattering length can 
be determined) and at somewhat higher energy where the J/ip-imdeon scattering amplitude 
may be determined. We use a model which takes into account quasi-elastic production from 
the deuteron (with production on one nucleon while the spectator nucleon recoils freely) as 
well as rescattering effects in the final state (proton-neutron rescattering and J/^/^-nucleon 
rescattering) . 

The motivation for the work in the first part of this paper (Sees. |n|-|lV|) was a proposal at 



JLab \Xl to measure the J/?/'- nucleon scattering length by the reaction 7* + — )■ J/ip +p + n. 
The proposed experiment would detect the outgoing proton and the decay products (e"'"e~) of 
the J/ip, for a kinematically complete measurement. The reason the J/'?/'-nucleon scattering 
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length is of interest is that several authors have argued that a nuclear bound state of the J 
may exist [31 S]. They propose that the force between & J /ip and a nucleon is purely gluonic in 
nature, and therefore is the analogue in QCD of the van der Waals force in electrodynamics, 
since the hadrons are color neutral objects. There is very little experimental data on elastic 
J/V'-nucleon scattering. There has only been one experimental measurement of it, at SLAC 
in 1977, where the J/'?/'-nucleon total cross-section was extracted by measuring production of 
J/ip's on nuclei and using an optical model for the re-scattering of the J/ip on the spectator 
nucleons [TB]. 

Measurement of the scattering length provides information on the bound states of the two 
particles involved in the scattering. In particular, for an attractive potential, if the scattering 
length is positive then there exists a bound state. Theoretical calculations of the J/'?/'-nucleon 
scattering length [3H6] predict a scattering length too small to produce a J/?/^- nucleon bound 
state. However, it would be large enough that there could exist a nucleus- J/?/^ bound state 
in a large enough nucleus. Since the scattering length is the (negative of) the zero-energy 
scattering amplitude, in order to measure this it is necessary for the two particles to scatter 
with small relative-momentum. In the case of 7* + — )■ J/ip +p + n at the energies which are 
kinematically allowed in the proposed JLab experiment, it isn't possible to have an on-mass- 
shell nucleon and J/ip scatter at small relative momentum. For an incident virtual photon 
of energy u = 9 GeV, and an outgoing J/'i/'-neutron pair with zero relative momentum, 
the minimum possible momentum of the neutron in the LAB frame (deuteron at rest) is 
~ 0.85 GeV; for u = 6.5 GeV and zero relative momentum of the J/V'-neutron pair, the 
minimum LAB momentum of the neutron is ~ 1 GeV. For zero relative momentum of the 
outgoing pair, the initial LAB momentum of an on-shell neutron in the deuteron (before the 
collision with the J/ip) must equal the final LAB momentum of the neutron. Therefore, the 
momentum of the neutron inside the deuteron would have to be 0.85 GeV (for u = 9 GeV). 
However, the deuteron wavefunction at that momentum is very small (essentially zero). 

So although the proposed experiment [1] may not be able to measure the J/T/^-nucleon 
scattering length, it might still be possible to measure the on-mass-shell J/'?/'-nucleon scattering 
amplitude, but at higher relative energies. The relative energy of the J/^/'-neutron pair would 
still be significantly smaller than in the only existing data (from the 1977 experiment at 
SLAC). Under certain kinematic conditions, the dominant contributions to the amplitude will 
come from p-n rescattering and/or J/ip — n rescattering after the J/tp is produced. If we fix the 
magnitude of the outgoing neutron's momentum at a moderately large value (here taken to be 
0.5 GeV) the contribution of the impulse diagram (where the J/tp is produced on the proton 
and the neutron recoils freely) will be negligible, since the impulse diagram is proportional to 
the value of the deuteron wavefunction at that momentum. This higher-energy rescattering 
is the subject of the second part of this paper (Sec. [V| ). 

This paper is organized as follows. Sec. [IT] reviews the general expression for the cross- 
section of electroproduction from a nucleus. In Sec. HI the kinematics for the case of zero 
and small relative momentum of the outgoing J/'?/'-neutron pair is discussed. In Sec. IV the 
calculation of the invariant amplitudes for 7* + /? — )■ J/ip +p + n are presented, including the 
one- loop diagrams corresponding to the p—n and J/ip-imdeon rescattering processes. In order 
to calculate the amplitude corresponding to the low-energy J/ip-neutron scattering, which 
involves the scattering length, model J/'j/'-neutron scattering wavefunctions and potentials 
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are used, and it is shown that the resuhing amphtude is insensitive to the model used. In 
addition, it is shown that for the kinematic conditions of the JLab experiment, the dominant 
amphtude is the impulse diagram, corresponding to J /ip production on the neutron with the 
proton recoiling freely, with no rescattering of any particles. This demonstrates that the 
measurement of the J/'?/'-nucleon scattering length is not feasible for the JLab experiment. 
Finally, in Sec. |V] calculations of the amplitude for 7* + D — )■ J/ip + p + n are presented 
under different kinematic conditions (not restricting the outgoing J/-?/;- neutron pair to small 
relative momentum). There it is shown that if the J/'j/'-neutron elastic scattering amplitude 
is somewhat larger than the value measured at SLAC at higher energy, it may be possible to 
extract this amplitude from the JLab experiment. 



II. ELECTROPRODUCTION FROM A NUCLEUS 

We consider here electron scattering from the deuteron with production of a vector meson, 
with the final state of the proton-neutron system being a continuum state. The formalism for 
the cross-section for electroproduction from a nucleus can be found in [9] . Here we summarize 
the relevant facts. We consider here the completely unpolarized electron (initial and final) 
cross-section. Then the cross-section can be written in terms of the amplitudes for 7* + — >■ 
p + n + V , i.e. vector meson production from virtual photons. In the LAB frame, with e' the 
final electron energy , Q' the final electron solid angle, fly the vector meson solid angle , and 
p*„ the final proton-neutron relative momentum in the p — n center-of-mass frame, the 8-fold 
differential cross-section has the form 

de'dn'dnyd^p* " (kinematic factors) x {vrR^i + vttR]J + vlR% + vtlR^i), (1) 
where the factors 



i?jr = 2Re(/|j;i(q)K)(/|J_i(q)K) (3) 



R^.t = -2Re(/|p*(q)K)((/|J+i(q)|.) - (/|J-i(q)K)) (4) 



i?f. = l(/|p(q)K)r, (5) 



are in terms of the matrix elements of the spherical vector components of the electromagnetic 
current operator J between the initial deuteron state \i) and final hadron {p + n + J/ip) state 
I/), vt, vtt, etc., are kinematic factors that only depend on the electron momenta. Rj^ is 
the sum of the squares of the amplitudes for '-f* + d^p + n + V for transversely polarized 
virtual photons, whereas R'j^ is an interference term between these two amplitudes. R^^ is 
the square of the amplitude for a longitudinally polarized virtual photon, while Rj^ is an 
interference term between the amplitudes for production from transverse and longitudinally 
polarized photons. The matrix element i^J/" is proportional to cos 4>, while R'j^ is proportional 
to cos 20, where is the angle between the plane including the initial and final electron 
momenta, and the plane including the 3-momentum transfer q and the J/ip momentum py. 
Thus if we integrate the cross-section over 0, the terms R^^ and R^^ drop out. Or, if we 
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assume helicity conservation (i.e. the helicity of the outgoing J/%1} is equal to the hehcity of 
the photon) then R^^ = R^^ = 0. Moreover, several theoretical models [TOl Hi] indicate 
that for small Q^, the amplitude for J/ip electroproduction from transverse virtual photons is 
much larger than the amplitude for production from longitudinally polarized virtual photons; 
for = (photoproduction) the production amplitude for longitudinal photon polarization 
is of course exactly zero. Therefore in what follows we will neglect R^^, and so the differential 
cross-section is simply given by multiplied by kinematic factors. Thus our task is to 

CcLlcU-lcttjG 

4. = K/l^+i(q)N)P + \{f\J-im)\' = |i^+r + \F-\' (6) 

where F± are the amplitudes for J/ip production from positive and negative helicity virtual 
photons. In the following we will calculate the amplitude for 'y* + d^p + n + Vhj evaluating 
Feynman diagrams corresponding to the various processes contributing to it. 

III. KINEMATICS FOR SMALL J/V'-NEUTRON RELATIVE MOMENTUM 

We will assume kinematics where the outgoing J/ip and neutron have a small relative 
momentum. Since the scattering length is the zero-energy limit of the scattering amplitude, 
in order to measure it the relative momentum of the J/ ip—n system must be small. An estimate 
of how small can be obtained by requiring only S-wave scattering, meaning the contribution 
of higher partial waves should be negligible. The classical relation between impact parameter 
and angular momentum yields an estimate for the maximum I that contributes. If the relative 
momentum of the J/tp — n pair is p* and the impact parameter is b, then the orbital angular 
momentum in the J/ip — n cm. frame is L ~ p*b = h^J2^T* where T* is the total kinetic 
energy in the 2-body cm. frame and /i is the reduced mass. The largest angular momentum 
partial wave which will be scattered is obtained by setting h equal to the range of the potential. 
With = l{l + 1) (we take h = 1), the condition for only S-wave scattering is that / -C 1, 
which implies = b'^2fiT* ^ 1. Taking the range of the interaction to be ~ 1 fm yields 
T* < 30MeV. 

Experimentally, perhaps the simplest quantity to measure is the total production cross- 
section, integrated over all available phase space, for a given incident photon energy. However, 
if we restrict the photon energy such that the maximum J/'i/'-neutron cm. kinetic energy T^^^ 
is small in of all the available phase space, then that means that the maximum proton- neutron 
relative energy will also be small everywhere in the available phase space: for a given value of 
the Mandlestam variable s for a system consisting of 3 particles, the total kinetic energy of 
any two of the particles (say 1 and 2) in their cm. frame satisfies 

T*2 < ^/s - nil - 1712 - m^. (7) 

The low-energy J/ip — n scattering amplitude is expected to be much smaller than the low- 
energy p — n scattering amplitude, and therefore the p — n rescattering would dominate over 
the J/ip — n rescattering, as contributions to the total production cross-section. Thus we 
need to restrict our considerations to a kinematic range where the p — n rescattering is at 
relatively high energy, while the J /ip — n rescattering is at very low energy, in order to have 
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FIG. 1: Momenta and angles in the overall cm. frame, for coplanar kinematics. 

the possibility that the J/ip — n rescattering makes a noticeable contribution to the differential 
cross-section. 

The ideal situation would be to have the final J /ip and neutron sitting at rest in the LAB, 
with the proton moving off at high velocity. Such a final state is kinematically allowed for other 
reactions, e.g. vr+d — )• rjpp, 7*6? — )■ r^pn, but it is not possible for the reaction 7*(i J /ip pn, 
for any real or virtual photon 4-momentum. 

Fig. |2] shows the minimum possible outgoing neutron LAB momentum and the correspond- 
ing proton LAB momentum vs. 9cm for fixed Ty^ = 30 MeV (see Fig. [ijfor the definition of 
6cm)', graphs are shown for photon LAB energy u = 9 GeV and for u = 6.5 GeV. One can 
see that the neutron's momentum is always greater than at least 0.6 GeV. Since the maxi- 
mum nucleon momentum in the deuteron is around 0.3 GeV (the deuteron momentum-space 
wavefunction is negligible for momenta larger than that) that means that in order for these 
final-state kinematics to occur, the neutron must have acquired its large momentum through 
a scattering event. In fact it will turn out that the dominant process corresponds to the 
impulse approximation wherein the J/ip is produced on the neutron itself, and the proton 
simply recoils freely. For the kinematics of interest here, rescattering processes (e.g. J/ip- 
neutron rescattering, J / ip-pioton rescattering, proton-neutron rescattering) make very small 
contributions to the total amplitude. 

We use the following notation throughout this paper: q = {u, q) is the virtual photon 
4-momentum in the LAB, with = —Q^ < 0; Pp is the outgoing proton LAB 3-momentum; 
Pn is the outgoing neutron LAB 3-momentum; pv is the J/ip LAB 3-momentum; and the 
same variables with cm superscripts denote their values in the overall (3-body) center-of-mass 
frame. 9p, 9n, and 9v denote the angle that the outgoing proton, neutron, and J/ip momenta, 
respectively, make with q, in the LAB frame. 
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FIG. 2: Minimum possible neutron LAB momentum (solid curve), and the corresponding proton 
LAB momentum (dashed curve), vs. 9cm, for Ty^ = 30 MeV and two values of photon LAB energy 
v. 




(c)^3a: J/tp-'O- rescattering diagram 



FIG. 3: Feynman diagrams for 7* + d — )• J/^ + p + n, for production on the proton. 



IV. INVARIANT SCATTERING AMPLITUDES 



The Feynman diagrams considered here are shown in Figs. [3] and |4j There are 3 dia- 
grams for production on the proton, and 3 similar diagrams where the J/ip is produced on 
the neutron. In all cases we are interested in kinematics where the J/ip and neutron have 
small relative momentum. The diagrams are covariant, and hence give Lorentz invariant 
amplitudes. In the diagrams, Ai'^^ is the Lorentz invariant amplitude for the quasi-2-body 
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{C')F3b- J/ip-p rescattering diagram 



FIG. 4: Feynman diagrams for 7* + d — )• J/ip + p + n, for production on the neutron. 

process 'j* + N ^ V + N (where is a nucleoli, and V stands for the J/ip)-, while A^^" 
is the Lorentz invariant amplitude for the elastic scattering process V + n — > V + n (with 
n meaning neutron), and Ai^^ and Ai"'^ are the same for elastic J/ip-pToton scattering and 
neutron-proton scattering, respectively. 



A. Impulse Diagrams 

Amphtudes Fia and Fib are the impulse diagrams, where the J/ip is produced on one of 
the nucleons and the other nucleon (the "spectator") recoils freely without interacting with 
the other particles. In Fia the vector meson is produced on the proton and the neutron is the 
spectator, while in Fib the production occurs on the neutron and the proton is the spectator. 
The invariant amplitudes in this case are 

Fla = M''''{Sia,tia) (8) 

Fib = M^''isib,tib)'^^. (9) 

Here Ai'^^ is the Lorentz invariant amplitude for the quasi- 2-body process 'y* + N ^ V + N 
(where is a nucleon, and V stands for the J /ip)-, To is the covariant vertex function for the 
virtual dissociation D — )■ p-|-n, and D{p) is the propagator denominator for the intermediate- 
state nucleon, D{j)) = —p'^ + rn? —ie. Si^, tia, Si^, and tib are the Mandlestam variables for the 
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2-body production process 'y* + N ^ V + N. Evaluated in the LAB frame, and neglecting 
any contributions to the deuteron vertex from antinucleons, the deuteron vertex function is 
related to the nonrelativistic deuteron wavefunction by [17] 

where in the LAB frame, k^e/ = P = — Pn for Fia, and krei = Pp = — n for Fib O^rei is 
the proton's momentum inside the deuteron, in the LAB frame, for both). In terms of the 
deuteron wavefunction, the amplitudes are thus: 



Sib, tlb) 

V'D(Pp)V2m(27r)3. (12) 

The amplitudes Ai"'^ used in calculations are to be taken from experimental data on J/ if) 
production on a single nucleon; we assume here that the amplitudes for production from a 
neutron is the same as for production from a proton. 

In the above expressions for the amplitudes Fia and Fib, spin labels have been suppressed. 
The initial virtual photon and the deuteron are in specific spin states, the final hadrons are in 
specific spin states, and there is a sum over the spin states of the intermediate-state virtual 
nucleon. For example, the amphtude Fib, including spin state specification, is explicitly: 

Fib = - V2m(27r)3 ^ A^^^(mi, A, m„, Xy) ^^'(pp, rm, m^), (13) 

nil 

where mi is the spin state of the intermediate-state neutron (i.e. the line with momentum n 
in the Feynman diagram), m„ and rrip are the spin states of the final neutron and proton, A is 
the photon polarization. Ay is the J/ip polarization, and M is the deuteron spin state. In what 
follows we will assume that the 2-body amplitudes for spin flip are negligible compared to the 
non-spin flip amplitudes, and so the amplitudes are diagonal in the nucleon spin, and also in 
the photon and J/ip spin. In that case we are able to calculate the spin-averaged squares of 
the various amplitudes Fia, F^a, etc., and the spin-averaged square of the total amplitude. 
We have included the contribution from the Z)-state in the deuteron wavefunction. For u = 9 
GeV the D-state was found to not make a significant contribution to the amplitudes, but for 
u = 6.5 GeV the D-state did contribute significantly, especially for the impulse diagram Fib. 
The deuteron wavefunction used was the Argonne f 18 wavefunction. 

All 2-body amplitudes Ai'^^ , Ai are related to the corresponding 2-body differential cross- 
section by 

^ = ^ \M\^ (14) 

dt 16irX{s,ml,ml) 

where the flux factor A is given in terms of the incident particle masses mi and m2 by 

A(s, ml, ml) = (s — m\ — m\Y — ^m\m\, (15) 



and s and t are the Mandelstam variables for the 2-body process. 
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1. Parameterization of the amplitudes M.^^ 
If the cross-section for J /ip production on a single nucleon is parametrized as 

with the parameters Ai, Bi dependent on energy (in principle), then the elementary produc- 
tion amplitude MC^ is given by 

M"^ = -I v/l67rAiA(s,-g2,m2)e5^i*, (17) 

where s and t are either s\a, t\a or s\h-, tib- 

The parameters Ai and Bi that are needed for the elementary J /if) production amplitude 
Ai'^'^ were estimated from the (scant) existing data on exclusive J/ip production on a nucleon. 
The only available data for the incident photon energy i/ ~ 10 GeV is from a photoproduction 
experiment at Cornell in 1975 [13]. For u in the range 9.3 to 10.4 GeV, they determined 
Ai = 1.1 ± 0.17 nb/GeV^ = (2.8 ± 0.43) x 10"^ GeV-^ and Bi = 1.31 ± 0.19 GeV-^. Those 
are the values used in this analysis. 



2. Fia and F^, 

Since Fia is proportional to \E'(p„), and as seen in Fig. [2]the outgoing neutron's momentum 
is always greater than ^ 0.6 GeV , the amplitude Fia will therefore be very small, since the 
deuteron wavef unction is negligible for those values of momentum. The amplitude Fib, on the 
other hand, is proportional to \I/(Pp); thus as seen in Fig. |2]for 9cm < 0.3 rad Fn, should be 
non-negligible since the proton momentum is less than 0.4 GeV over that range of 9cm- 



B. One-loop diagrams 

The covariant expression for a general one-loop diagram (see Fig. [s] ) is 

r d'n Tjn') M^'^'M 
J i{2TiY D{n') D{n)D{k) ^ ' 

where n', n, k are the internal 4-momenta indicated in the figure, M. stands for either 
A^P", A1^", or (elastic scattering amplitude for proton-neutron, V-neutron, or V-proton 

scattering, respectively), and D{j)) = — nip + ie, etc., are propagator denominators. Spin 
labels have been suppressed in Eq. ( jlsj ); in particular, there is an implicit sum over the spin 
states of the intermediate-state particles (the lines labelled n, k, and n'). There are 4 diagrams 
total for a given set of outgoing proton, neutron and J/ip momenta. Taking p2 = p„, pi = Pp 
(so that the internal line n is the neutron, and n' and k are the proton) gives one diagram 
(p-n rescattering diagram). The other 3 are: p2 = Pn, Pi = Pv, where the internal line n is 
the neutron, k is the J/ip, and n' and ps are the proton (V-n rescattering diagram); p2 = Pp, 
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FIG. 5: General one- loop diagram, n and p2 are the same particle (either neutron or proton). The 
line n can be either a proton or a neutron. 

Pi = Pn, where the internal line n is the proton, and n' and k are the neutron (another p-n 
rescattering diagram); and P2 = Pp, Pi = Py, where the internal line n is the proton, k is the 
J /ipi and n' and p3 are the neutron (V-p rescattering diagram). 

All of the one-loop diagrams can be evaluated in the same manner. We follow closely the 
method of j^. We first integrate over by taking the positive-energy pole at vP = ujn — 
where a;„ = \/w? + n^, coming from the propagator denominator D{n). This corresponds to 
neglecting antinucleon components of the deuteron wavefunction. We also use the relation 



Eq. (10) (and we work in the rest frame of the deuteron). This gives 



-^(n) . (19) 



(27r)3 J 2ujn D{k) 

Note that in this expression, the internal nucleon line n is now on-mass-shell, since n° = Un- 
Thus in the time-ordered diagram of Fig. [Sj only the lines n' and k can be off-shell; all the 
rest are on-shell. 

The above expression for the amplitude F can be separated into two terms, one term in 
which the line k is on-mass-shell and one in which k is off-mass-shell, by using the relation 

^ -i7i6{x) + V- (20) 



X + ie X 



for the remaining propagator in Eq. (19) (with V representing the principal value) and by 
choosing the coordinate system shown in Fig. [6] where the Uz axis is along the direction 
of the vector pi2 = pi + p2. The delta function term gives the on-mass-shell part of the 
amplitude, and the principal value term gives the off-mass-shell part. This separation is useful 
since the elementary amplitudes Ai"^^ , Ai can be determined (at least in principle) directly 
from experimental data on the relevant 2-body scattering processes only when all 4 particles 
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FIG. 6: Coordinate system used. For the on-shell amplitude, 9 is fixed for a given n = |n| and pi, 

P2- 

involved (2 initial and 2 final) are on-mass- shell. This is not the case if one of the particles 
involved is off-mass-shell. But if for reasonable choices of the off-mass-shell amplitudes the 
off-mass-shell part is small compared to the on-shell part, then the off-mass-shell part will 
not play an important role. The result for the on-shell and off-shell parts of the one-loop 
amplitude is [7]: 

F = + F"^^ (21) 



where 



and 



-tTT 



v/2m(27r)3 2|pi2| 



I d(f) dnn'^{n)M'^^M (22) 

Jo J\n-\ 



Here 



F"f^ = —= , / d<\) / dnn^(n) V / c/cos0— — -. (23) 

v/2m(27r)3 2|pi2| Jo Jo J /i2(n) + cos^ ^ ' 

2n |pi2| 

with n = |n|. Si2 = (^1+^2)^ is the Mandelstam s- variable for the elastic scattering of particles 
1 and 2, rrik is the mass of the real particle which the line k represents, £"12 = Ei + E2, and 6 
is as shown in Fig. [6j Note that /12 is independent of 9 and 0. 

In the amplitude the limits of integration |n_| and n+ are the solutions of 

/i2(n±)2 = 1, (25) 

which are 

n± = ^|pi2|±^£i2 (26) 



where P2, E2 are the momentum and energy of outgoing particle 2 in the cm. frame of 
particles 1 and 2, and particle 2 is the same particle as the internal line with momentum n. 
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The range of n given by |n_| < n < ri+ is the range of n for which it is kinematically possible 
for the hne k to be on-mass-shell (given that n, pi, and p2 are on-shell). In the value of 
cos 9 is fixed at 

cos 6 = -fi2{n) = . (27) 

2|pi2F 

The amplitudes M^^ , M in F"" are evaluated, for a given n and 0, at this value of cos 6. The 
amplitude A4 in is fully on-shell, i.e. all 4 particle lines n, k, pi, and p2 are on- mass-shell. 
The amplitude Ai"^^ has only one particle off-shell (n'), but given that the magnitude of n is 
small (due to the deuteron wavefunction), n' is almost on-shell: n'^ = Ma — Un — Mfi — m ~ m, 
and so (n')2 = m2 + C(g). 

In the amplitude F"^^ , k is never on-mass-shell; the principal value imposes this, since for 
k to be on-mass-shell, cos 6' must equal —fuin), which never occurs in the principal value. 
Thus the amplitudes Ai'^^ , Ai that enter into have either one particle off-mass-shell 
(for Ai) or two particles off-shell (for Ai"'^). In our calculations the amplitude F"-^-'^ was 
much smaller than F"", for any of the Feynman diagrams in Fig. [s] and |4] (except for Fsa), 
for the kinematics considered here, and for reasonable expressions for the off-shell amplitudes 
M^^, M. 



C. General features of the one-loop diagrams 



Inspection of Eq. (22) shows that a necessary condition for the on-shell part of a given one- 



loop diagram to be non- negligible is that the corresponding \n_\ must be small enough so that 



the range of integration in Eq. (22) includes the momenta where the deuteron wavefunction 



is significant. Since n_ depends on the final-state kinematics (Eq. (26)), graphs of n_ as a 
function of one of the final state variables reveal regions where the on-shell amplitude can be 
significant, and regions where it will be negligible. 

For the diagram of most interest, F^a, particles 1 and 2 are the J/ip and neutron, respec- 
tively. Fig. [T] shows n± vs. 9cm for Ty^ = 30 MeV, for i/ = 9 GeV, and Fig. 7(b) shows the 



same for z/ = 6.5 GeV, for the diagram F^a- At both of these photon energies, n_ is greater 
than ~ 0.6 GeV for all 6cm, and so the on-shell amplitude will be negligible since the deuteron 
wavefunction is negligible for that momentum. 

For the amplitudes F2a, F2b and Fsf,, the corresponding n_ graphs are shown in Fig. [9] for 
z/ = 9 GeV and z/ = 6.5 GeV, for Ty^ = 30 MeV. Note that at a given value of 6cm and 
Ty„, the neutron LAB momentum can range from a minimum to a maximum allowed value 
(with these two values both corresponding to p„ and pv pointing in the same direction in the 
LAB, with 6*' = in Fig. [s]), and the values of for F2a, F2b and F^b depend on the neutron 
LAB momentum in addition to 6cm and Ty^. For our calculations, we have fixed the neutron 
LAB momentum for a given 6cm (and Ty^ = 30 MeV) at its minimum value; we denote this 
value by Pn,min- One can see from these graphs that for z^ = 9 GeV, there are intervals of the 
variable 6cm for which the on-shell parts of F2a, F2b, and F^b should be non-negligible, since 
n_ < 0.05 GeV there (note that for the diagram F^b-, the J/i/^-proton rescattering occurs at 
relatively high energy, if the J / ^/'-neutron relative energy is small; so F^b is not directly related 
to the J / ip-micXeon scattering length). For v = 6.5 GeV, n_ is larger than ~ 0.4 GeV, and so 
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FIG. 7: n± vs. 6cm for the amplitude F^a (J/V'-neutron rescattering) , for Ty.^ = 30 MeV, and 
I' = 9.0 GeV and u = 6.5 GeV. Upper curve is n_(_, lower curve is n_. 




FIG. 8: LAB frame momenta. For a given 6cm and Ty^, the neutron LAB momentum can range 
between a minimum and maximum value, both of which correspond to 6' = in the figure. PnV = 

PV + Pn- 

these on-shell amplitudes should be small. This is born out by the exact calculations, where 
the one-loop on-shell amplitudes for u = 6.5 GeV are in general much smaller than those for 
u = 9 GeV. 



D. Parameters used in the elementary 2-body amplitudes 



For the calculation of the amplitudes, using Eq. (22) and Eq. (23), the elementary am- 
plitudes Ai'^^ and Ai are phenomenological amplitudes obtained from existing experimental 
data, related to the 2-body cross-section by Eq. (14). We take the individual 2-body differen- 
tial cross-sections to be of the form 
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(a)z/ = 9 GeV 
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FIG. 9: n_ vs. 6 

T^^ = 30 MeV, pr, 



for proton-neutron rescattering and J/V'-proton rescattering amplitudes, for 

Pn,min- 



where the A's and 5's can depend on energy. Thus we have 



M 



(29) 



relating to A and B. 

The values of A and B depend on the relative momentum (or energy) of the rescattering 
pair. Table [T] lists the values of the momentum p of the neutron in the proton's rest frame 
(for the p — n subsystem) and the momentum p of the J/ip in the proton's rest frame (for the 
V — p subsystem), for the case of Ty^ = 0. Note that these quantities are independent of 9cm 
(easily shown in the overall cm. frame). Also given in Table |l] is the total kinetic energy of 
the pair in the cm. frame of that pair, T*2- Given the values of the momentum p, we can 
determine the parameters that enter into the elementary amplitudes A^^", Ai^^. 



TABLE I: Parameters used in elementary scattering amplitudes 
1/ (GeV) Subsystem p (GeV) (GeV) B (GeV-^) atot (mb) A (GeV-^) 
9 n-p 2.25 0.64 5.7 - 6.2 43 - 46 260 

9 V-p 7.4 1.02 1.31 3.5 1.61 

6.5 n-p 0.86 0.16 6.9 35 160 

6.5 V-p 2.84 0.247 1.31 3.5 1.61 



1. p — n scattering parameters 
For z/ = 9 GeV, the existing data [S] for p — n scattering at incident momentum p = 2.25 



GeV gives B 



pn 



5.7 to 6.2 GeV- 



The value of Apn can be obtained from the total 



p — n cross-section by using the optical theorem and neglecting the real part of the scattering 
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amplitude: 

^l*=o=^a,, = A,„ (30) 

The measured value of atot given in the table is then used to calculate Apn. For u = 6.5 GeV, 
the existing data for p — n scattering at momentum p = 0.86 GeV gives Bpn = 6.9 GeV~'^ [8] 
and Apn = 160 Gel/"^ [I2]. 



2. J ftp — p scattering parameters 

There is very little data on elastic J/ip-pioton scattering from which to determine the 
parameters Avp and Bvp that are needed for the J/'?/'-proton rescattering amplitude M^^. 
For the present analysis, we have assumed that the t-slope for elastic J/'?/'-nucleon scattering 
is equal to the t-slope for the process 7* + — )■ J/ip + N, and so we've taken Byp = 
B^v = 1-31 ± 0.19 GeV^^. We can obtain Ayp from the total J/V'-nucleon cross-section, 
using the optical theorem; however, there has only been one measurement of (jjl^ ^ [T5] . 



In an experiment in 1977 at SLAG [16] J/ip photoproduction was measured on beryllium 
and tantalum targets, and the total J / tp-micXeon cross-section was extracted by using an 
optical model for the rescattering of the produced J/V' from the other nucleons in the nucleus. 
The value they obtained was afl^ ^ = 3.5 ± 0.8 mb, which gives via the optical theorem 
Ayp = 1.61 ± 0.4 GeV~'^. In that paper, however, they also note that the measured J/ip- 
photoproduction cross-section together with vector meson dominance arguments would give 
a J/ip-imdeon total cross-section of ~ 1 mb. So we can assume the value of the J/'j/'-nucleon 
total cross-section to be not very well known. In addition, the photon energy in the SLAG 
experiment was 20 GeV, and so assuming forward production of the J/ip, then the energy 
of the J/ip in the LAB frame would also be ~ 20 GeV, giving a kinetic energy in the LAB 
of ~ 17 GeV. This is significantly larger than the kinetic energy of the J/ip in the proton 
rest frame considered here, where for z/ = 9 GeV it is 4.94 GeV and for u = 6.5 GeV it is 
1.1 GeV. This introduces more uncertainty in the value of Ayp to be used. In a theoretical 
calculation of the J/ip-imdeon scattering length yields a value for the total J / ip-imdeoia cross- 
section at threshold of 7 mb, and it is argued that the total cross-section should decrease as 
the energy is increased from threshold. Thus at the energy of the J/V'-proton rescattering 
here, the value of Ayp may be larger than the value measured in the experiment at SLAG. 
For the purpose of calculating the amplitude F^b, however, we will use the value measured at 
SLAG. 



3. Subthreshold J/ip production 
The threshold photon energy for production on a single nucleon at rest is 

l^thresh = my^ (31) 

2m 
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while for production on the deuteron it is 

J^thresh = my H — (32) 

For Q2 = 0.5 GeV, these are 8.47 GeV and 5.78 GeV, respectively. 

For u = 6.5 GeV, which is below threshold for J/ip production on a single nucleon at rest, 
we assume that the production mechanism is the same as for production on a free nucleon. 
The Fermi motion of the nucleon in the deuteron is what allows the production to occur, i.e. 
if the nucleon is moving towards the photon with a large enough momentum then the value 
of Si = {q + p)"^, where p is the 4-momentum of the nucleon in the deuteron, will be above 
the threshold value. In the calculation of the amplitudes for i/ = 6.5 GeV this condition was 
imposed on the internal nucleon momentum in the integrals involved. 



E. Calculated On-shell and Off-shell amplitudes 

Using the parameters in Table |T], the on-shell and off-shell parts of the amplitudes were 
calculated. The squares of the individual amplitudes F2a, -^26, and F^h are shown in Figs. 



10 and 11 shown in the graphs is a curve which includes only the (square of the) on-shell 
part of the amplitude, and also a curve which is the square of the total amplitudes including 
both the on-shell and off-shell parts. For the off-shell parts, the same parametrizations of the 
elementary amplitudes Ai^^ and Ai were used as for the on-shell parts. As stated previously, 
the off-shell parts are very small compared to the on-shell parts, which means that knowledge 
of the exact forms of the off-shell elementary amplitudes A4'^^ and Ai are not needed. 

Since the J/ip-imdeon scattering length is expected to be small (much smaller than e.g. 
the proton- neutron scattering length), the J/'i/'-neutron rescattering diagram F^a should be a 
small contribution to the total amplitude. This is born out in the next subsection, where 
is calculated using a model potential and wavefunction, for a value of the scattering length on 
the order of that predicted by theoretical models. 



F. J/^-neutron Rescattering Diagrams and the Scattering Length 



Diagram F^a (Fig. 3(c)[ ) is the J/'i/'-neutron rescattering diagram. This is the diagram 
where the J/ip and neutron scatter from each other with small relative momentum; hence this 
amplitude will involve the scattering length for the J/'i/'-neutron interaction. 

For this diagram, we have pi2 = Pvn = Pv + Pn in Eqs. (22) - (27), and 9 is the angle 
between py„ and n (see Fig. |6|. For T* = 0, where T* = T* + Ty {* indicates the cm. frame 
of the V — n system), we have = |n_| (see Eq. (26)), and so F^^ = 0. For small T* (below 
e.g. 100 MeV) 



0, where T* = T* + [* 
, ...... - ^see Eq. (|26[), and so F"^" 

will be small because n_ > 0.45 GeV for the possible JLab kinematics 
see Fig. [7|. Thus the main contribution to F^a is from F^j,^ , for which the intermediate-state 



Tpon 
^3a 



J/ip is always off-mass-shell. However, because of the propagator denominator fuin) + cos 6*, 
where 

1 



fuin) + cose 



2\pvn\n 



(33) 
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contributions to -^3/'^ from values of n for which k is far off-mass-shell will be small. To obtain 
estimates of F^l^ , we will therefore evaluate it using on-mass-shell values of Ai"'^ and Ai^^. 

The relation between the invariant amplitude A^^" and the scattering amplitude f{k,6) 
is [H] 

M = -8n^f{k,e) (34) 

for the on-energy-shell amplitudes. The half-off-energy shell amplitudes are related in the 
same way. 

Our normalization conventions are the following. The scattering length a is given by 

\imf{k,e) = -a. (35) 

fc— >o 

The off-energy-shell amplitude is given by 

r^(ki,k2) = -(2vr)V (k2|V^|<^) = -(27r)V (<^|V^|ki) = -(2vr)V (ki|y|<V (36) 
where V is the potential, /i is the reduced mass, ki is the initial relative momentum (in terms 



of n, k in Fig. 3(c)), k2 is the final relative momentum (in terms of py, p„ in Fig. 3(c)), 



and \E'k2 is the exact scattering wavefunction for asymptotic relative momentum k2. Since 
this off-energy-shell scattering amplitude depends on the scattering wavefunction "^^J and 



\F2„f (GeV-2) IFibf (GeV-2) 




6cm (rad) 



(c) 



FIG. 10: Squares of amplitudes for u = 9.0 GeV, T^„ = 30 MeV, 

Pn — Pn,min- The solid curve is 

the total (on-shell plus off-shell parts), while the dashed curve is only including the on-shell part of 
the amplitude. 
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the potential V, both of which are unknown for J/t/'-nucleon elastic scattering, we will resort 
to models in order to estimate the amplitude. 

We normalize our 5"- wave scattering wavefunction ^^k, and define the radial wavefunction 
u{r), by: 



Mr) 



1 



iS{k) 



u{r) 



(37) 



7(2^ 

where 5{k) is the 5- wave phase shift. In order to calculate the matrix element, we will specify a 
model zero-energy wavefunction UQ{r) which determines the potential V{r) via the Schrodinger 
equation, and use that potential to solve for the wavefunction for A; 7^ (the subcript zero on 
Uq indicates it is for A; = 0). 

We assume the J / ip-rmclcon potential is of finite range, and so is zero for r larger than 
some distance R. The phase-shift 5{k) satisfies the following well-known properties as k — > 0: 

1. for a repulsive potential, or an attractive potential that doesn't admit a bound state: 
5 — )■ —ak as A; ^ 0; 

2. for an attractive potential which admits a single bound state: 5— >7r — a/cas/c— >0 
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FIG. 11: Squares of amplitudes for u = 6.5 GeV , T*„ = 30 MeV, 

Pn — Pn,min- The solid curve is 

the total (on-shell plus off-shell parts), while the dashed curve is only including the on-shell part of 
the amplitude. 
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The zero-energy wavefunction \l/o(^') outside the range of the potential is then 



r — a 



3 



for both cases, while the zero-energy radial wavefunction Mq"* differs by a minus sign for the 
two cases; this is purely due to including the factor e'^^^'^^ in the definition of \l'fc in Eq. (37). 



Our normalization conventions give a > for either a repulsive potential or an attractive 
potential with a bound state, and a < for an attractive potential that doesn't admit a 
bound state. In all cases a is the intercept on the r-axis of u^uti''^)- 

Theoretical calculations [HE] give values of \a\ ~ 0.3 /m, with effective range — 2.0 fm. 
It is thought that the potential is attractive, but too weak to support a bound state. Below, 
calculations of F^a are made for both cases of an attractive potential: a > (bound state) 
and a < (no bound state). 



G. Positive scattering length a 

For the case of a positive scattering length and attractive potential (which possesses a 
bound state), we have 'UQ"*(r) = — (r — a) . We choose the simplest zero-energy wavefunction 
Uo{r) consistent with the standard continuity requirements on u{r) at r = and r = R 
required by the Schrodinger equation. This yields: 



One further requirement on Mq" is that Mq" have no zeros on the interval [0, R] (besides at 
r = 0). This ensures that the corresponding potential V{r) is non-singular, since for k = 0, 
V{r) = 2^^- The zeros of Mq" are at r = and 



r 3 1 R 

and the right-hand-side must lie outside the range [0, 1]. This requires either 

R<a (41) 

or 

R>3a (42) 

The potential is 

11 r - R 

Vir) = — - - - f43l 

^ ' 2fir R^{3- ^) ~3Rr + r^ ^ ' 

and one can see that for R> 3a the potential is repulsive. Therefore we require R < a. 
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1. Model wavefunction for k ^ 

Given this model potential we can proceed to calculate the off-shell scattering amplitude 
Eq. (36) and the amplitude F^f'^ once we calculate the wavefunction u for non-zero k for a 



given a and R. Taking a = 0.3 fm, R = 0.1 fm, and solving the Schrodinger equation for 
k < 100 MeV numerically, we find that the wavefunction \E'k('") and the off-shell amplitude 
f{p, k) depend very weakly on k over the entire momentum range from k = to k = 100 
MeV. 

Since we are only interested in the J/^/'-neutron relative momentum up to around 100 
MeV, it is legitimate to approximate the off-energy-shell amplitude f^"'{p,k) ~ /^"(p, 0) for 
the range of k we are interested in. For our model wavefunction Mq" we can evaluate /^"(p, 0) 
analytically: 

/^"(p,0) = / dVe-P-f/(r)vl/o(r) 

V (2^) (44) 
6a / sin pR'^ 



■1 + 



p2i?2 V pR 

The momentum p appearing in f{p,0) is the relative momentum of the J/V'-neutron pair in 
their center-of-mass frame, before they scatter in diagram F^a] it is thus the magnitude of n 
(or -k) in the outgoing V — n center-of-mass frame, and so we must boost n to that frame. 



H. Negative scattering length a 



have UQ^'i^r 



1 



(45) 



If the potential is attractive but too weak to support a bound state, then a < and we 
r — a . In this case the zero-energy wavefunction is 

-V + f-V^-f 

RJ \RV \R^J 

The requirement that u have no zeros on [0, R] imposes no restriction on a and R in this 
case. Theoretical calculations give a around —0.3 fm, and effective range ~ 2.0 fm jHE]. 
Using these values with our model wavefunction and potential implies i? = 1.3 fm. Again 
there's very little variation with k of the wavefunctions and off-shell amplitudes for k from 
to 0.1 GeV, so to calculate F^a in this case we will again approximate /^"(p, k) ~ /^"(p, 0). 
This yields 

y„, 6a sin pi?" 



/^"(P,0) 



pR 



(46) 



Results 



The electroproduction differential cross-section in the LAB frame, shown in Figs. 12 and 
13 , is given by 



voVt E' 



pI 



dE'dVL'dppdVLpdVLn S{2TifMd E 8(27r)5 Ep \Ev pi - EnPn ■ P 



v\ 



(47) 
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FIG. 12: Squares of amplitudes, and electroproduction differential cross-section Eq. (47), for v = 
9 GeV, r*„ = 30 MeV, 

Pn — Pn,min- The solid curves in (b) and (c) includes all amplitudes, while 
the dashed curves (which are not distinguishable from the solid curve) omit i<3a for both positive 



and negative scattering length a, for the model potential of Eq. (43). 



where F = Fi„ + Fib + + Fat + Fg^ 4 
a virtual photon, and 

yieF^F 



^^0 



12 



and 



Vt 



E' = E 



Fsfe is the total amplitude for J/ip production from 

(48) 
(49) 



2 q2 



+ 



4FF' cos2(^72), 
^0 ' 



12 GeV 



V. 



(50) 

In the above, E is the initial electron energy (taken to be 12 GeV), E' is the final electron 
energy, and Q' is the scattering angle of the electron relative to the initial electron momentum 
(all quantities in the LAB frame). 

The results are shown in Figs. 12 and 13 For u = 9 GeV, the squares of the individual 
amplitudes Fu, F2b, F^b, and F^a are also shown in Fig. 12(a), along with the square of 



the total amplitude including all diagrams. As can be seen in that graph, the amplitude 
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FIG. 13: Electroproduction differential cross-section Eq. (47) for v = 6.5 GeV , Ty^ = 30 MeV, 
Pn = Pn,min, and a = ±0.3 fm. In (c), the dashed curves (not distinguishable from solid curve) are 
omitting for both positive and negative a, for the model potential of Eq. (43). 



F^a makes a negligible contribution to the total amplitude. There are intervals of 6cm where 
the amplitudes Fn,, F2b, and ^3;, individually dominate the total amplitude. However, by 
comparing Fig. 12(a) with Fig. 12(c), which shows the electroproduction differential cross- 
section on a linear scale, over the range of 6cm for which the cross-section is non- negligible (for 
< 6cm < 0.2 rad) the cross-section is due exclusively to the impulse diagram Fib. The very 
small "bump" visible in Fig. 12(c) at 6cm — 0.3 rad is due to the proton-neutron rescattering 
amplitude F2b- 

For u = 6.5 GeV, the difference between the cross-section including F^a and omitting ^3^ 



is visible in the logarithmic-scale graphs (Figs. 13(a) and 13(b)), but not in the linear-scale 



graph. Fig. 13(c) 

F^a was also calculated using 3 other potentials: a square-well potential yielding a = 0.3 fm 
and R = 0.1 fm, and also the potential of Eq. (43) but with a = 0.3 fm, R = 0.29 fm, and 
a square-well potential yielding a = 0.3 fm but with R = 0.29 fm. There was a negligible 
difference between the values of ^3^ calculated with these potentials. 
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J. Conclusion 

It does not appear to be possible to measure the J/^/'-nucleon scattering length via pro- 
duction on the deuteron, under the kinematic conditions available at JLab. For small values 
of the relative momentum of the outgoing J/i/j-neutron pair, the initial momentum of the 
neutron inside the deuteron that is required for on-mass-shell rescattering of the J/?/'- neutron 
pair is larger than ~ 0.6 GeV (see Fig. [?]), where the deuteron wavefunction is negligible. 
The off-mass-shell part of the rescattering amplitude was calculated using model J/'?/'-nucleon 
potentials and was found to make a negligible contribution to the total amplitude. The vast 
majority of J/ip production events, for Ty^ < 0.03 GeV, will be at small values of 9cm, 
where the impulse diagram Fn, dominates, and therefore information on J/^-nucleon elastic 
scattering at small relative energy cannot be obtained. 



V. INTERMEDIATE ENERGY J/^j PRODUCTION ON THE DEUTERON 

It may be possible to extract the J/ip-nudeon elastic scattering amplitude from the 7* + 
D — !■ J/i{j+p+n experiment, at higher relative energy of the J/-?/'-nucleon pair, under different 
kinematic conditions for the final-state particles than was considered in the previous sections of 
this chapter. Under certain kinematic conditions, the dominant contributions to the amplitude 
will come from rescattering diagrams (p-n rescattering and J/ip — n rescattering). If we fix 
the magnitude of the outgoing neutron's momentum at a moderately large value (here taken 
to be 0.5 GeV) the contribution of the impulse diagram will be negligible, since the impulse 
diagram is proportional to the value of the deuteron wavefunction at that momentum (see 
Fig. [4] for the impulse and rescattering diagrams). For the analysis presented here, we: 

• use coplanar kinematics 

• fix the magnitude of the outgoing neutron momentum at Pn = 0.5 GeV 

• fix the 4-momentum-transfer-squared t = (g — py)^ at a particular value 

• plot amplitudes or differential cross-sections vs. ^„ (the angle that the outgoing neutron 
momentum p„ makes with direction of the incoming photon momentum) for fixed Pn 
and t 



see Fig. 14). For some range of z/ and t, these graphs will display peaks due to p — n and 



J /if) — n on-mass-shell rescattering. For z/ = 10 GeV , the peak due to J /ip — n rescattering 



is evident (see Fig. 17), but for i/ = 9 GeV it is not evident (see Fig. 18). This analysis is 
similar to what has been done in [7] for the reaction '-y + D-^T^ + N + N. 

The kinematics here are very different than in the previous sections. There it was the 
relative energy of the J/'?/'-neutron system which was kept fixed, at a small value, while the 
parameter which was varied was the production angle of the J /il) in the overall center-of-mass 
system. 
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FIG. 14: (a) "plus" kinematics and (b) "minus" kinematics. For "plus", Pp is always on the opposite 
side of the photon momentum q direction as the neutron momentum. For "minus" , py is always on 
the opposite side of the photon momentum q direction as the neutron momentum. 



A. Intermediate-energy J/ip production 



We are interested here in kinematics available at JLab after the 12 GeV upgrade. The 
maximum (virtual) photon energy is then around 11 GeV. Here we evaluate the amplitude 
for virtual photon 4-momentum q = {u, q) with u = 10 GeV, = —q^ = 0.5 GeV, keeping 
Pn = 0.5 GeV and t fixed. The results presented here are for t = —2 GeV^; calculations were 
done for larger values of \t\, with similar results (although the total amplitude decreases with 
increasing \t\). We consider the same set of diagrams as before, shown in Figs. [3]and|4j The 
impulse diagrams, Fia and Fn,, are negligible for these kinematics. 

Note that for a given t, pn, and 6n, there are 2 sets of allowed values of the proton and 
J/ip momentum {pp,Pv}; we've called the two sets the "plus" set and the "minus" set. If we 
define x and z axes as in Fig. [14} with the a;-component of the neutron momentum always 

and the 



positive, then the "plus" kinematics is as shown in Fig. 14(a 



■'minus" kinematics is 



as shown in Fig. |14(b) For the "plus" kinematics, Pp^ is negative for all 6'„ (while pvx takes 



both positive and negative values over the range of 9n), while for the "minus" kinematics, pvx 



is negative for all 9r, 
Figs. 15 and Fig. 



(while Ppx takes both positive and negative values over the range of On)- 
16 show graphs of |n_| vs. On for the 3 different pairs of outgoing 
varies greatly with On, the value of the amplitude of the 



n_ 



particles. Since the value of 
corresponding rescattering diagram varies greatly also, and has a prominent peak at the value 
of On for which |n_| = 0. 



Calculation of amplitudes 



For the calculation of the amplitudes, the elementary 2-body amplitudes Ai'^^ and Ai 
are taken to be of the diffractive form Aea^* with parameters determined from existing 
experimental data. For the J/'?/'-nucleon rescattering diagrams, the only available data is 
from the experiment at SLAG [12] discussed in Sec. IV D They determined the total J/ip- 
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FIG. 15: |n_| vs. for photon energy u = 10 GeV, and t = —2 GeV^. (a) and (b) are for the 
"minus" kinematics, (c) and (d) are for the "plus" kinematics. 
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(a) J/V'-p rescattering (b) J/'^-p rescattering 



FIG. 16: |n_| vs. 0„ for J/tp-P rescattering, for photon energy u = 10 GeV, and t = —2 GeV"^. (a) 
is for the "minus" kinematics, (b) is for the "plus" kinematics. 

nucleoli cross-section to be cyfl^ ^ — 3.5 ± 0.8 mh, which gives via the optical theorem 
Avn = 1-61 ± 0.4 GeV^'^. The energy of the J/ip in this experiment was ~ 20 GeV in 
the Lab frame (nucleon at rest). However, for our kinematics the rescattering of the J/ip on 
the nucleon takes place at an energy in the outgoing neutron's rest frame of from 6 to 10 GeV, 
which is significantly smaller than in the SLAG experiment; thus the value of Avn at our 
energy may be significantly different. Since the entire reason for measuring the cross-section 
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for this process is to extract the J/V'-nucleon scattering amphtude in an energy region where 
it has not been measured before, we've used several different values of the parameter Avn in 
the calculations, from 1 times the SLAG value up to 10 times the SLAG value. Since the total 
cross-section atot for J/'j/'-nucleon scattering goes like y/ Avn (Eq- (30)), this corresponds to 
a range of atot (which is what was actually measured in the SLAG experiment) of from 1 to 
~ 3 times the SLAG value. In theoretical calculation of the J/'i/'-nucleon elastic scattering 
cross-section at threshold yielded 7 mb, which is twice the value measured at the higher energy 
at SLAG. 

The full calculation of the amplitudes must of course include the off-shell parts. If we 
use the same parametrization of the elementary amplitudes as in the on-shell part, then the 
off-shell parts were found to be very small compared to the on-shell parts. 

The amplitudes and F^t, where the J/ip is produced on the neutron and then rescat- 
tering (of the neutron or J/ip, respectively) occurs on the proton, are much smaller than 



and F^a, and do not exhibit the well-defined peaks that F2a and F^a do. Fig. 17 shows the 



8-fold electroproduction differential cross-section, Eq. (47), versus 9n- In that figure, the neg- 
ative values of 6'„ are for the "minus" kinematics, and the positive values are for the "plus" 
kinematics. Graphs are shown for 3 different values of the J/?/'- neutron elastic scattering pa- 
rameter Avn'- Avn = 1-6 GeV~'^ (which is the value determined in the experiment at SLAG), 
Avn = 8.0 GeV~^, and Avn = 16 GeV~^. It is seen that only if Avn is of the order of 10 times 
as large as the previously measured value is there a noticeable peak due to the J/V'- neutron 
rescattering, for the "plus" kinematics. The p — n rescattering peak is much larger than, and 
close enough to, the J/-?/'-neutron rescattering peak that it obscures the J/ip peak. For the 
"minus" kinematics, the same statement holds; in addition, however, the size of the p — n 
rescattering peak varies (by ~ 40%) as the value of Ayn is varied. Note also that the position 
of each of the peaks is simply given by the value of 6'„ where the value of the corresponding 
|n_| is zero (see Fig. 15). 

It's important to note that the peak due to the J/^/'-neutron rescattering isn't observable 



at lower energies. Fig. 18 shows the square of the total amplitude for photon energy of 



u = 9 GeV and t = —3, for Avn = 10 x 1.6 GeV On this graph the peak due to p — n 
rescattering is visible, but there's no visible peak due to J/?/)- neutron rescattering. 



B. Conclusion 



We have shown here the possibility of measuring the elastic J/ ?/^-nucleon scattering ampli- 
tude for energies significantly smaller than the energy of the only existing data. If the total 
J /ip-nncleon cross-section o'H^^ at these energies is of the order of 2 — 3 times the previ- 
ously measured value, then the differential cross-section as a function of On should exhibit 
well-defined peaks corresponding to on-mass-shell p — n and J /if) — n rescattering, for virtual 
photon energy of = 10 GeV and 4-momentum-transfer-squared t = {q — pvY = ~2 GeV^. 
However, at lower photon energy (below 9 GeV) the J/ip — n rescattering peak would not be 
distinguishable. As it is expected [1] that oHf^ should increase as the energy decreases, it 
is not impossible that the lower energy cross-section could be larger than the measured value 
by a factor of ~ 2. 
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FIG. 17: Electroproduction differential cross-section vs. including all diagrams, for photon energy 
1/ = 10 GeV, Q2 = 0.5 GeV"^ and t = -2 GeV'^, for 3 values of Avn- Solid curve: Avn = 1-6 GeF"^. 
Dashed curve: = 8.0 GeV~^. Dotted curve: Avn = 16 GeV~^. In (a), the large peaks at 

9n — 45° and 9n ~ —85° are due to p — n rescattering, while the small peaks (or bumps) at On — 80° 
and 9n — —70° are due to J/V'-neutron rescattering. (b) shows detail of left half of (a). 



IFI^GeV-^) 




FIG. 18: Amplitude squared vs. 9n, including all diagrams, for Avn = 16 GeV^^, for photon energy 
u = 9 GeV and t = —3 GeV^. The peak or bump on either side is due to p — n rescattering, while 
the peak due to J/V'-nucleon rescattering is not visible. 



VI. CONCLUSION 



This paper is concerned with determining information regarding J/T/'-nucleon elastic scat- 
tering. As discussed in Sec. IV J, it does not appear possible to measure the J / ip-nudeon 



scattering length via production on the deuteron under the kinematic conditions available 



at JLab. As discussed in Sec. VB, however, it will be possible to extract the J/V'-nucleon 



scattering amplitude at higher relative momentum if (tH^^ is large enough. 
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